Cooperativity and the origins of rapid, single-exponential kinetics in protein folding 
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The folding of naturally occurring, single domain proteins is usually well-described 
as a simple, single exponential process lacking significant trapped states. Here we fur- 
ther explore the hypothesis that the smooth energy landscape this implies, and the 
rapid kinetics it engenders, arises due to the extraordinary thermodynamic coopera- 
tivity of protein folding. Studying Miyazawa-Jernigan lattice polymers we find that, 
even under conditions where the folding energy landscape is relatively optimized (de- 
signed sequences folding at their temperature of maximum folding rate), the folding 
of protein-like heteropolymers is accelerated when their thermodynamic cooperativity 
enhanced by enhancing the non-additivity of their energy potentials. At lower temper- 
atures, where kinetic traps presumably play a more significant role in defining folding 
rates, we observe still greater cooperativity-induced acceleration. Consistent with 
these observations, we find that the folding kinetics of our computational models more 
closely approximate single-exponential behavior as their cooperativity approaches op- 
timal levels. These observations suggest that the rapid folding of naturally occurring 
proteins is, at least in part, consequences of their remarkably cooperative folding. 

PACS numbers: 87.15. Cc; 91.45.Ty 

Keywords: Monte Carlo simulation; optimal folding temperature, native topology, contact or- 
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INTRODUCTION 

Diffusive processes, such as protein folding, progress 
more rapidly across a smooth energy landscape than 
on rough landscapes dominated by myriads of energetic 
barriers (Bryngelson 1995). Exhaustive lattice and off- 
lattice heteropolymer simulations support this simple ar- 
gument; rapid folding is observed only among those rel- 
atively rare sequences that lack well-populated, trapped 
intermediates (Mirny 1996; Chan and Dill 1997; Chan 
and Dill 1998; Cieplak 1998). Similarly supporting this 
argument, the rapid folding of naturally occurring, single 
domain proteins is almost always associated with energy 
landscapes that are smooth relative to fcsT, as indicated 
by the observation of simple, single exponential kinetics 
even at the lowest temperatures that can be investigated 
experimentally (Gillespie and Plaxco 2000; 2004). 

While the relationship between smooth energy land- 
scapes and rapid folding is well established (Bryngel- 
son 1995; Mirny 1996; Chan and DiU 1997; Chan and 
Dill 1998; Cieplak 1998), the perhaps more subtle ques- 
tion of the molecular origins of this smoothness has seen 
relatively less attention. One of the few attempts to 
address this issue directly is the work of Shakhnovich 
and co-workers, who have noted that the landscapes of 
the vast majority of lattice polymer sequences are very 
rough, and that smooth lattice polymer landscapes are 
comparatively uncommon (Miller 2002). Based on this 
observation they speculate that polypeptides exhibiting 
smooth landscapes are similarly uncommon and that the 
rapid folding of naturally occurring proteins arises as a 



consequence of evolutionary selections aimed at ensuring 
that this critical property is achieved (e.g. refs Abkevich 
1996; Gutin 1998a; Mirny 1998). Unfortunately, how- 
ever, experimental studies have not supported this pro- 
posal; the folding of de novo designed proteins -proteins 
produced without any regard to the roughness of their en- 
ergy landscapes- is generally rapid (Gillespie 2003; Zhu 
2003; Scalley-Kim and Baker 2004; Khulman and Baker 
2004), suggesting that smooth landscapes are a rela- 
tively common property of thermodynamically foldable 
polypcptides(Gillespie 2003; Gillespie and Plaxco, 2004). 
Thus the ultimate origins of the smooth folding energy 
landscapes observed for naturally occurring proteins re- 
mains unclear. 

It has been speculated that the discrepancy between 
the rapid, trap-free folding observed for most small pro- 
teins and the generally slow, trap-dominated folding of 
lattice and off-lattice models (Jewett 2003) may arise be- 
cause, in contrast to the folding of most computational 
models, protein folding is thermodynamically coopera- 
tive (Kaya and Chan 2003; Kaya and Chan 2003a; Chan 
2004; Kaya 2005). Thermodynamic cooperativity, usu- 
ally identified by the calorimetric criterion of Privalov 
(see e.g. Makhatadze and Privalov 1995), is defined as 
a bimodal conformational population peaked around the 
native and unfolded state enthalpies and practically zero 
at intermediate enthalpies. The word cooperativity is 
used to describe this effect because it is widely thought 
to arise from nonadditive (i.e. cooperative) interactions 
akin to those first proposed as an explanation of the co- 
operative oxygen binding of hemoglobin (Wyman and 
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Allen 1959). The speculation that similar cooperativity 
(arising from similarly non-additive interactions) might 
account for the smooth landscapes almost universally ob- 
served for the folding of small proteins is based on the 
hypothesis that such will cooperativity accelerate folding 
more by destabilizing misfolded, trapped states than it 
decelerates folding by destabilizing productive interme- 
diates (Jewctt 2003; Eastwood and Wolynes 2001). 

It has not proven possible to modulate the thermo- 
dynamic cooperativity of proteins in the laboratory, and 
thus it has not yet proven feasible to explore the relation- 
ship between cooperativity and folding rates in vitro. It 
is possible, however, to modulate the cooperativity of 
computational protein models and thus simulation ex- 
periments provide a useful means of testing the hypoth- 
esized link. Despite extensive simulations-based studies 
of the role of thermodynamic cooperativity in generating 
linear chevron behavior, topology-dependent rates and 
other signatures of two-state folding kinetics, however, 
the extent to which cooperativity acts to accelerate or de- 
celerate folding has not been explicitly investigated in the 
prior literature (Fan 2001; Shimizu and Chan 2002; Kaya 
and Chan 2003; Jewett 2003; Kaya and Chan 2003a; 
Kaya and Chan 2003b; Cieplak 2004). Here we explore 
this question via simulations of protein-like Miyazawa- 
Jernigan lattice polymers to which increasing degrees of 
cooperativity have been introduced. 

RESULTS 

Test sysytems 

We have explored the folding behavior of 48-residue 
lattice polymer sequences encoding two topologically dis- 
tinct native structures. The two structures, named 
topologies 1 and 2, exhibit relative contact orders (Plaxco 
1998a) of 13% and 26% respectively, placing them among 
the least and most "complex" topologies that a maxi- 
mally compact, 48-monomer structure can adopt (Figure 
1; table 1). Using the design algorithm of Shakhnovich 
and Gutin (Shakhnovich and Gutin 1993) we obtained 
three, rapidly folding sequences encoding each of these 
distinct topologies. The sequences were designed to 
adopt the target topologies using the Myiazawa-Jernigan 
(MJ) energy parameterization, 

JV 

E{W.}, in}) = E ^'^^ - (1) 

i>j 

where {fl} is the set of bead coordinates that define a 
conformation, {at} represents an amino acid sequence, 
and ai stands for the chemical identity of bead i. The 
contact function A is 1 if beads i and j form a contact 
(that is not a covalent linkage) and is otherwise. The 



energy parameters Eij are taken from the 20 x 20 MJ 
matrix, derived from the distribution of contacts of na- 
tive proteins (Miyazawa and Jernigan 1985). In order to 
minimize the possiblity of energy-related kinetic effects, 
sequences were selected that exhibit similar native state 
energies, Enat (table 2). 



Enhancing cooperativity 

While it is clear that a remarkable feature like ther- 
modynamic cooperativity demands highly unusual ener- 
getics, the nature of the interactions underlying the ther- 
modynamic cooperativity of protein folding has not yet 
been determined (Takada 1999; Kaya and Chan 2000; 
Fernandez 2002; Shimizu and Chan 2002; Kaya 2005). It 
is thought, however, that the thermodynamic coopera- 
tivity of protein folding arises due to multi-body effects 
leading to nonadditive interactions (i.e. in which the for- 
mation of one "bond" makes the formation of subsequent 
"bonds" more favorable (Wyman and Allen 1959)), a hy- 
pothesis that is supported by extensive simulation stud- 
ies (Eastwood and Wolynes 2001; Fernandez 2002; Kaya 
and Chan 2003b; Ejtehadi 2004). In order to capture the 
nonadditivity that presumably underlies thermodynamic 
cooperativity we have employed a modified version of the 
MJ potential that entails many-body, nonpairwise inter- 
actions in the form of a nonlinear relationship between 
the energy of a conformation E and the fraction of native 
contacts, Q, it forms: 

where 5 is a parameter measuring the magnitude of 
the induced cooperativity and E is given by Equation [1] 
Of note, Q is the fraction of native contacts; when we 
employ a similar function that counts the total number 
of contacts we find that MJ polymers are not thermo- 
dynamically stable if S rises even trivially above unity. 
Note too that when S ~ 1 we recover Equation 1 and 
the traditional, noncooperative MJ potential. As S in- 
creases, however. Equation 2 promotes a steeper decrease 
of the protein's energy with increasing number of na- 
tive contacts when folding nears completion. It does 
so by penalizing the energy of conformations in which 
only a small fraction of the total set of native contacts 
is formed {Q <~ 0.7; data not shown) and by "repay- 
ing" that penalty after that threshold is crossed. As ex- 
pected, when S is increased as per equation [2] the ther- 
modynamic cooperativity of our model systems is sig- 
nificantly enhanced; both the enthalpy distribution and 
conformational distribution of our model polymers be- 
come strongly biomodal (i.e. see e.g. Makhatadze and 
Privalov 1995) as S is increased from to Sopt (Figure 2) 
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FIG. 1: In this study we have employed sequences folding into two topologies (numbered 1 and 2), which are among the least 
and most complex topologies attainable for a maximally compact 48-residue lattice polymer. (Bottom) The bands along the 
main diagonal in the contact maps (Saitoh 1993) of the two topologies bands indicating contacts between one amino acid and its 
four successors represent the structural equivalent of a-helices. /3-sheets are represented as thick bands parallel or anti-paxallel 
to the diagonal. 



The kinetic consequences of cooperativity at the optimal 
folding temperature 

In order to address the extent to which enhanced co- 
operativity enhances folding rates wc initially performed 
simulations at T^pt, the temperature at which the folding 
rate of each sequence reaches its maximum value in the 
absence of added cooperativity. Under these conditions 
we find that folding invariably accelerates as cooperativ- 
ity, S, increases towards its optimal level, Sopt' at their 
Sopt the folding rates of all six of the sequences wc have 
characterized increase by 1.14 to 3.2-fold relative to the 



rates observed at 5 = 1 (Figure 3, left). Above Sopt fold- 
ing rates decrease, presumably because further increases 
in cooperativity significantly destabilize the folding tran- 
sition state. 

While the folding time of every test system wc have 
investigated exhibits a similar functional dependence on 
S, we observe quantitatively different behavior among 
differing structures and sequences. For example, under 
these conditions Sopt ranges from 1.2 to 1.3 for the three 
sequences that fold into topology 1 and the associated 
rate increases (over the non-cooperative, S' = 1 bench- 
mark) range from 2.0 to 3.2-fold. In contrast, over a 
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FIG. 2: Enhancing the nonadditivity (i.e. cooperativity -sec Wymann and Allen, 1959) of the MJ energy potential enhances 
the thermodynamic cooperativity (Makhatadze and Privalov 1995) of MJ lattice polymer folding. As illustrated here using two 
representative sequences at their folding transition temperatures, T/, the energy distributions (top row) and conformational 
distributions (bottom row) of our model polymers become more strongly biomodal (i.e. more thermodynamically cooperative - 
see e.g. Makhatadze and Privalov 1995) as S is increased from 1 to Sopt- In particular, the population of fully native molecules 
(bottom row, Q = 1.0) is enhanced significantly at higher values of S. 



similar range of Sopt we observe only 1.14 to 1.8-fold ac- 
celerations for sequences folding into the more complex 
topology 2. 



The kinetic consequences of cooperativity at lower 
temperatures 

If cooperativity accelerates folding by smoothing the 
energy landscape, it might be expected to have a larger 



effect at temperatures below Top* , where landscape 
roughness plays an increasingly important role in defining 
folding kinetics (Gutin 1996; Gutin 1998a) (under these 
conditions there is a higher probability for the chain to 
get trapped in metastable states). This is perhaps all the 
more true for sequences which, like those we have em- 
ployed, were designed using an algorithm ensures their 
energy landscapes are relatively smooth (Shakhnovich 
and Gutin 1993). Consistent with this expectation we 
find that the effects of induced cooperativity are indeed 
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FIG. 3: The folding of MJ lattice polymers is accelerated when the cooperativity of the MJ energy potential is enhanced. 
This effect holds for all six sequences we have investigated (adopting both low -top row- and high -bottom row- contact order 
structures). And while this acceleration is readily apparent even at the temperature of optimal folding (Topt; left column), it is 
significantly more pronounced at lower temperatures (T = Topt - 0.03; right column), where kinetic traps might be expected to 
play a more significant role in defining kinetics. Above some optimal level of cooperativity {Sopt), however, folding decelerates 
with increasing cooperativity. This presumably occurs because further increases in S destabilize native elements in the folding 
transition state, slowing folding more than the destabilization of trapped states accelerates it. 



much more striking at temperatures below Tgpt : at T 
= (Topt - 0.03) we observe up to 6-fold increases in fold- 
ing rate (Figure 3, right). The coopcrativity-induced in- 
crease in folding rates is, in fact, so great at T = (Topt 
- 0.03) that under these conditions (at S = Sopt ) all 
six sequences fold more rapidly than they do at Topt 



(for any S). This presumably occurs because Topt rep- 
resents a compromise between the driving force behind 
folding (native state stability), which increases as T de- 
creases, and the kinetic consequences of trapped states, 
which tend to slow folding when T is reduced. Because of 
this, lowering the temperature below Topt and raising S 
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FIG. 4: At the still lower temperature of T = Topt - 0.05 
(corresponding to a ~ 16° C tcmpcratnrc drop for a protein 
with a Topt ^ 40° C e.g. the similarly sized fynSfiS doniain51) 
the accelerating effects of enhanced cooperativity are even 
more pronounced. Under these conditions the folding rate 
of a sequences encoding topology 1 (top) increase by 5.7 to 
6. 9- fold as S increases from unity to Sopt- This is more than 
twice the maximum enhancement observed at Topt and 2 to 
40% greater than the enhancement observed at T = (Topt 
- 0.03). For a single sequence adopting the more complex 
topology 2 (bottom) we observe a 4.4-fold acceleration at this 
lower temperature, which is 80% greater than that observed 
at Topt and 28% greater than that observed at T = (Topt - 
0.03). 



to Sopt may accelerate folding by increasing the driving 
force (Gutin 1996; Gutin 1998b; Shakhnovich 1994) with- 
out concomitantly stabilizing trapped states that would 
otherwise slow the process. Perhaps consistent with this 
argument, it has previously been noted that rapidly fold- 
ing lattice polymers tend to be those associated with a 
large gap between the energy of the native state and that 
of all other maximally compact states (Sali et al., 1994), 
an effect that might also arise due to a relationship be- 
tween folding kinetics and thermodynamic cooperativity 
(Chan 2004). To date, however, it has not yet been shown 
whether such an energy gap leads to thermodynamic co- 
operativity. 

At the still lower temperature of 0.05 energy-units be- 
low Topt (which corresponds to a ~ 16°C temperature 
drop for a protein with a Topt of ^ 40°C; sec, e.g., ref 
Plaxco 1998) folding is significantly slowed and thus the 
exhaustive simulation of all six of our test sequences be- 
comes difficult (particularly for sequences folding to the 
more complex topology 2). Nevertheless, we have investi- 
gated all three topology 1 sequences and a representative 
topology 2 sequence (sequence 2) under these conditions 
and find that the cooperativity-induced acceleration con- 
tinues to increase as the temperature is reduced. For 
example, at this much lower temperature the three se- 
quences encoding topology 1 fold 5.7 to 6.9 times more 
rapidly at Sopt than at S = 1 (Figure 4, left). These 
accelerations are more than twice those observed at Topt 
and 40 to 200% greater than those observed at T = (Topt 
- 0.03). Similarly, the folding of the sequence adopting 
the more complex topology 2 increases 4.4-fold at this 
temperature (Figure 4, right), which is 1.8 times the ac- 
celeration observed at Tgpt and 28% greater than that 
observed at T = (Topt - 0-03). These accelerations are 
so significant that, even at this very low temperature 
(which, if for no other reason, should slow folding due to 
the increase in the barrier height relative to ksT), the 
folding of the cooperativity-optimized heteropolymers is 
significantly faster than that observed at Topt in the ab- 
sence of added cooperativity. We note, however, that 
although it is typical of this class of lattice models (see 
e.g. Faisca 2002a, Faisca 2005) the dispersion of fold- 
ing times we observe is smaller than that observed for 
single domain proteins (which span a six order of magni- 
tude range). The chains we have investigated, however, 
arc homogeneous in both length and stability, which per- 
haps accounts for the relatively limited dispersion in their 
folding rates 

Cooperativity enhances rates by smoothing the energy 
landscape 

If thermodynamic cooperativity accelerates folding by 

reducing the stability of trapped states (i.e. smoothing 
the energy landscape), it might be expected to push ki- 
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netics from heterogeneous, non-exponential behavior to- 
wards simpler, single-exponential behavior (Bryngelson 
and Wolynes 1987; Onuchic 1997; Socci 1998; Nymeyer 
1998). Consistent with this, Chan and Kaya have shown 
that significant thermodynamic cooperativity must be 
built into Go -type lattice polymers in order to recover 
the linear chevron behavior observed for real proteins 
(Kaya and Chan 2003a). Using two sequences (adopting 
topologies 1 and 2 respectively) we have investigated this 
question in more detail (Figures 5). We find that, while 
the folding of these sequences are generally rather well 
approximated as single exponential (because they were 
designed to exhibit smooth energy landscapes), without 
exception the kinetics of our models nevertheless move 
closer to single-exponential behavior as S increases to 
its optimal value. For example, in the absence of added 
cooperativity we observe correlation coefficients (for the 
proportional relationship between log(fraction unfolded) 
versus the number of MC steps) ranging from = 0.950 
to 0.976. At Sopt, in contrast, these correlation coeffi- 
cients invariably increase, ranging from 0.988 to 0.998. 
Thus modest levels of added cooperativity appear to in- 
variably improve the single-exponential behavior of even 
well designed sequences, suggesting that the added coop- 
erativity is accelerating folding by smoothing the folding 
energy landscape. 



Cooperativity may enhance the topology- dependence of 
folding rates 

We (Jewett 2003), and others (e.g. Kaya and Chan 

2003a; Ejtehadi 2004), have previously noted that modifi- 
cation of the Go potential to enhance its cooperativity en- 
hances the extent to which Go polymer folding rates cor- 
relate with measiires of native state topology. Given that 
we have investigated only two topologies in this study 
(albeit they are among the most and least topologically 
complex of the maximally compact 48-mers), we cannot 
similarly establish whether increasing cooperativity also 
enhances the topology-dependence of MJ polymer folding 
rates. Nevertheless, our results are consistent with this 
speculation: whereas, at Topt, the three sequences that 
adopt the more complex (higher CO) topology 2 struc- 
ture fold, on average, 6.2-times more slowly than those 
adopting to topology 1, this discrepancy almost doubles 
to 11.9-fold when S rises to Sopt (Figure 3). Whether 
this effect truly reflects a cooperativity-induced enhance- 
ment of the already significant topology-dependence of 
MJ polymer folding rates (Faisca and Ball 2002a), how- 
ever, will require the simulation of significantly more 
topologies than were sampled here. 



DISCUSSION 

It is known that the folding of Go polymers, which are 
characterized by a native-centric energy potential and 
thus exhibit vc;ry smooth energy landscapes, slows as 
their cooperativity is increased (Jewett 2003; Kaya and 
Chan 2003a). Such increase in the folding time results 
from the fact that an increase in cooperativity leads, in 
such modified Go -type models, to a destabilization of the 
transition state relative to the original Go model which, 
by considering only native interactions, is explicitly en- 
ergetically biased towards the native state. Here we have 
shown that, in contrast, the folding rates of energetically 
more complex -and thus perhaps more realistic- MJ mod- 
els, where native as well as non-native interactions con- 
tribute to protein energetics, is increased when their ther- 
modynamic cooperativity is increased and that this ef- 
fect is enhanced as the temperature drops and landscape 
roughness plays a larger role in defining kinetics. The 
observed acceleration presumably arises because modest 
cooperativity increases folding rates more by destabiliz- 
ing misfolded, trapped states than it decelerates the pro- 
cess by destabilizing the folding transition state (East- 
wood and Wolynes 2001). In support of this claim, the 
folding kinetics of all six of the sequences we have investi- 
gated approach single-exponential behavior more closely 
as their cooperativity is increased and, as Chan and co- 
workers have noted (Kaya and Chan 2003a), increases 
in thermodynamic cooperativity tend to enhance linear 
chevron behavior. 

While our rescaling of the energy leads to a smoother 
energy landscape we should note that this does not neces- 
sarily imply that the potential we have employed becomes 
more Go-like as S is increased. Go and Go-like mod- 
els exhibit smooth energy landscapes because they are 
native-centric, i.e., they neglect the effect of non- native 
interactions and only consider the contribution of native 
interactions to the protein's energy. The cooperativity 
term we have introduced, in contrast, enhances the sta- 
bility of native-like conformations only when folding is 
near completion (i.e., for Q >0.7). Moreover, this effect 
is mild at the modest levels of cooperativity employed 
here {S <1.5) . For smaller values of Q our rescaling pro- 
cedure destabilizes the energy of any conformation with 
a certain fraction Q of native contacts irrespective if that 
conformation is en-route to folding or not. This means 
that the stability of kinetic traps is not necessarily and 
selectively reduced relative to that of native-like confor- 
mations. Indeed, would this be the case and, contrary 
to finding an optimal value for S, one would expect to 
observe a monatonic increase in folding rate with increas- 
ing S. The lack of such a relationship implies that the 
mechanism by which landscape smoothness is achieved in 
our model is differs from the mechanisms underlying the 
smoothness of Go models and that, in turn, the observed 
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accelerations are not a trivial outcome of the manner in 
which cooperativity has been encoded. 



The extent to which cooperativity accelerates folding 

appears to depend on both sequence-specific and topo- 
logical effects. Indeed, it appears that increasing cooper- 
ativity may increase the spread in folding rates between 
simple and more complex topologies (Kaya and Chan 
2003a; Jewett 2003; Ejtehadi 2004; Kaya 2005) in a man- 
ner reminiscent of the topology dependence of protein 
folding rates (Plaxco 1998a). This topology dependence 
presumably arises because, while the added cooperativity 
destabilizes energetically trapped, misfolded states (i.e. 
reduces so-called "energetic frustration"), it docs not af- 
fect the entropic barriers that arise from the polymer 
properties of the chain, such as connectivity (Makarov 
and Plaxco 2003) and excluded volume effects, and quirks 
of the native topology, such as lack of symmetry (Nel- 
son 1997) (together sometimes termed "topological frus- 
tration" (Shea 1999)). Such a linkage between folding 
rates, topology and cooperativity is consistent with the 
topomer search model (Makarov and Plaxco 2003), which 
postulates that two-state folding kinetics are defined by 
the diffusive search for unfolded conformations sharing 
the native topology (i.e. are in the "native topomer"). 
If folding is thermodynamically cooperative, only the na- 
tive topomer can form sufHcient native-like interactions 
to ensure productive folding and thus the entropic cost of 
finding the native topomer will dominate relative folding 
rates. Similarly, while lattice polymers cannot fold via a 
strictly topomer-search process (a coarse lattice polymer 
cannot be in the native topology without actually being 
in the native state), the simulations described here illus- 
trate the more general observation that, as cooperativity 
is increased, global properties such as native-state topol- 
ogy will play an increasingly important role in defining 
the folding barrier (Eastwood and Wolynes 2001; Plotkin 
1997). 



Irrespective of the origins of the effect, previous 
studies clearly demonstrate that increasing thermody- 
namic cooperativity significantly enhances the topology- 
dependence of lattice polymer folding rates (Kaya and 
Chan 2003a; Jewett 2003; Ejtehadi 2004). Here we have 
shown that the addition of cooperativity also acceler- 
ates the folding of energetically complex lattice polymers 
and enhances the single-exponential nature of their ki- 
netics. Taken together, these observations suggest that 
the observed topological dependence of two-state protein 
folding rates may be a consequence of the cooperativity 
necessary to ensure the smooth energy landscapes upon 
which rapid, single-exponential kinetics can occur. 



METHODS 

The lattice polymer model 

We consider a simple three-dimensional lattice model 
of a protein molecule with chain length L = 48. Individ- 
ual amino acids are represented in this model by beads of 
uniform size occupying lattice vertices. These monomer 
units are connected into a polypeptide chain using bonds 
of uniform (unitary) length corresponding to the lattice 
spacing. 

In order to mimic folding to the native state we use a 
standard Monte Carlo (MC) algorithm together with the 
kink-jump move set (see e.g. Landau and Binder 2000) 
and with the energy defined as in Equation 2. Accord- 
ingly, random local displacements of one or two beads are 
repeatedly accepted or rejected in accordance with the 
standard Metropolis rule (Metropolis 1953). All MC sim- 
ulations were started from randomly generated unfolded 
conformations and their folding dynamics were traced by 
following the evolution of the fraction of native contacts, 
Q = n/N, where N is the total number of native contacts 
and n is the number of native contacts formed at a given 
MC step. Each simulation was run until the polymer 
reached the native fold, i.e. when Q = I. The number 
of MC steps required to fold to the native state was de- 
fined as the first passage time and the folding time was 
computed as the mean first passage time (MFPT) of 100 
simulations. 

Structure selection and sequence design 

We selected two model topologies from a pool of 500 
maximally compact conformations (MCC) found using 
simulations of homopolymer relaxation. The latter are 
MC simulations where a polymeric chain composed by 
beads of the same chemical species is launched, at tem- 
perature T = 0.7, from a randomly generated conforma- 
tion and relaxes to the minimum energy conformation. 
For self-attractive homopolymers of length L = 48 on 
a three-dimensional cubic lattice, the most stable con- 
formation is a MCC cuboid with 57 contacts. We have 
computed a histogram of the contact order frequencies 
of a large sample of these MCCs (data not shown) and 
identified the two extremes at CO values of 13% and 
26%. These topologies, numbered here 1 and 2 respec- 
tively, were selected for our studies (Figure 1 and table 
!)• 

Ensembles sequences were designed to fold into topolo- 
gies 1 and 2 (table 2). The goal of the design process was 
to identify sequences that fold rapidly and efficiently into 
a pre-selected conformation named target structure. We 
followed Shakhnovich and Gutin (Shakhnovich and Gutin 
1993) and Shakhnovich (Shakhnovich 1994) in attempt- 
ing this by seeking the sequence with the lowest possible 
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TABLE L Geometric traits for native topologies 1 and 2. Contact order, CO, number of long-range contacts Qd, we define 
as long-range (LR) a contact between any two contacting beads i and j whose sequence separation is at least 12 units, i.e. 
\i — j\ > 12 (Gromiha and Selvaraj 2001; Makarov and Plaxco 2003). Long-range order, LRO, which is equivalent to Qd 
normalized by chain length (Gromiha and Selvaraj 2001) 

Topology # CO Qd LRO 

1 0.13 0.33 0.40 

2 0.26 0.74 0.88 



TABLE IL Designed sequences. The folding rate is measured at the optimal folding temperature Topt- 



# (Topology #) 


Sequence 




Topt 


T/ 


-logio(i) 


1 (1) 


FRTRPLNHDFYNYKIWEPFKPADFPKAWDRMLDHVWDSMASWGHQHCS 


-25.85 


0.32 


0.32 


-6.91 ± 0.05 


2 (1) 


CDLPPFTYRHHGNDFWKNYEMIKHWDLWRDMFRAFWSDPVKASPHQAS 


-25.92 


0.32 


0.32 


-6.42 ± 0.03 


3 (1) 


FRTPWVSHQFYAYKLMEHFKWGDFCRNMDKWIDSLPDRWNPAPHDHAS 


-26.09 


0.32 


0.32 


-6.18 ± 0.04 


1 (2) 


KDKIHFRMNYGYPAWDAQSVKDLTCPRDWHFPHMRDPSHNWELAFFWS 


-25.87 


0.34 


0.33 


-7.47 ± 0.05 


2 (2) 


ENDVTMDMDPSPCLFRIHNLPRAHSFDRFGWHQFDKYHYKWKWAWAPS 


-26.15 


0.31 


0.34 


-7.30 ± 0.04 


3 (2) 


EHDAQLDFDWSRWTWHGRNSYHAPAMYRWPVHDMDKPNPKFKIFFLCS 


-26.24 


0.32 


0.33 


-7.44 ± 0.06 



energy in the target state, as given by Equation 1. To this 
end the target's coordinates were quenched and the en- 
ergy of Equation 1 annealed with respect to the sequence 
variables a. This amounts to a simulated annealing in 
sequence space: starting from some randomly generated 
sequence, transitions between different sequences (which 
are generated by randomly permuting pairs of beads) are 
successively attempted along with a suitable annealing 
schedule. 



Optimal folding temperature 

The existence of optimal folding temperature (Topt), 
corresponding to the temperature at which the folding 
rate is maximized, has extensively reported in lattice 
studies of protein folding (Gutin 1996; Gutin 1998b; 
Shakhnovich 1994; Faisca and Ball 2002; Cieplak 1999; 
Faisca and Ball 2002a; Faisca 2005). In order to deter- 
mine Topt we have computed the folding times (in the 
absence of added cooperativity) over a broad tempera- 
ture range for all six of the sequences we have investi- 
gated. Figure 6 illustrates the dependence of the loga- 
rithmic folding rates, logio{l/MFPT), on the simulation 
temperature for sequences 1 and 2 folding to topologies 
1 and 2 respectively. These sequences exhibit Topt of 
0.32 and 0.31 respectively. Below and above these opti- 
mal temperatures their folding slows significantly, with 
the stronger temperature dependence being observed at 
lower temperatures. 

For comparison with prior lattice polymer folding stud- 
ies we have determined the folding transition tempera- 
tures, T/, of the six sequences we have employed. Tj 
-typically denoted as the melting temperature or T^ in 
the experimental literature) is the temperature at which 



denatured states and the native state are equally popu- 
lated at equilibrium. In the context of a lattice model it 
can be defined as temperature at which the average value 
< Q > of the fraction of native contacts is equal to 0.5 
(Abkevich 1995). In order to determine Tf we averaged 
Q, after collapse to the native state, over MC simulations 
lasting at least 20 times longer than the average folding 
time computed at Topt,- Similarly, in order to compute 
population histograms at T as a measure of thermody- 
namic cooperativity (Figure 2) data were averaged after 
collapse to the native state over MC simulations least 20 
times longer than the average folding time at Topt- 
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FIG. 5; As the cooperativity increases to its optimal value (Sopt) the folding of MJ lattice polymers more closely approximates 
a single-exponential relaxation. Shown are results for two representative sequences (sequence 1 of topology 1 and sequence 
2 of topology 2) over a range of temperatures. While the folding of these sequences is generally well approximated as single 
exponential (because they were designed to fold efficiently) , a distinct improvement in the single-exponential behavior of folding 



13 



-6.8 




0.26 0.28 0.3 0.32 0.34 0.36 

temperature (T) 

FIG. 6: The folding kinetics of MJ lattice polymers are 
strongly temperature dependent. Shown is the dependence of 
the logarithmic folding rate, logio{l/MFPT), on the simula- 
tion temperature T, for sequence 1 of topology 1 and sequence 
2 of topology 2, under conditions of no induced cooperativity 
{S = 1). The optimal folding temperature is the temperature 
at which folding is most rapid. 



